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We propose a 5-dimensional definition for the physical 4iI>-Yang-MiUs theory. The fifth 
^ . dimension corresponds to the Monte-Carlo time of numerical simulations of QCD4. The 
■ 5-dimensional theory is a well-defined topological quantum field theory that can be renor- 
malized at any given finite order of perturbation theory. The relation to non-perturbative 
physics is obtained by expressing the theory on a lattice, a la Wilson. The new fields that 
must be introduced in the context of a topological Yang-Mills theory have a simple lattice 
expression. We present a 5-dimensional critical limit for physical correlation functions and 
for dynamical auto-correlations, which allows new Monte-Carlo algorithm based on the 
time-step in lattice units given hj e = Qq ^^/^^ in pure gluodynamics. The gauge- fixing in 
five dimensions is such that no Gribov ambiguity occurs. The weight is strictly positive, 
because all ghost fields have parabolic propagators and yield trivial determinants. We 
indicate how our 5-dimensional description of the Yang-Mills theory may be extended to 
fermions. 



1. Introduction 

The 4-dimensional Yang-Mills theory seems to suffer from logical contradictions. In 
the continuum formulation, one has a gauge-fixed BRST invariant path integral, but one 
has the famous Gribov ambiguity for large gauge field configurations [1][2]. One often 
discards this problem, since the idea of defining the theory as a path integral of a gauge field 
can only be seriously advocated in perturbation theory or for semi-classical approximations. 
In the lattice formulation, which is by construction valid non-pertubatively, one chooses as 
variables the gauge group elements, but one has yet another contradiction. The way the 
continuum theory is approached is unclear, and if one tries to do a local gauge-fixing, the 
partition function vanishes [3] . This question is also often discarded, since for computing 
gauge-invariant quantities, one can factorize the volume of the gauge group that is finite 
on the lattice. 

It is however frustrating not to have continuum and lattice formulations which would 
separately define both gauge-invariant and gauge non-invariant sectors, with a BRST sym- 
metry controlling the correspondence of the theory to a physical sector, and with a natural 
limit from the lattice to the continuum formulation. 

In this article, we will show that, in order to reconcile the continuum and lattice 
approaches, it is useful to define the theory in a 5-dimensional space, such that the 4- 
dimensional physical theory lives in a slice of this extended space. The theory that we 
will consider passes the tests that are obviously needed: in the continuous formulation it 
is perturbatively renormalizable by power counting without loosing its physical character, 
(due to Ward identities), and the gauge- fixing is no longer subject to the Gribov ambiguity; 
in the lattice formulation, all fields that one considers in the continuum formulation take 
their place, and one now obtains a formulation with a consistent gauge fixing. 

The fifth-dimension will be the stochastic time that Parisi and Wu proposed long 
time ago for stochastically quantizing the Yang-Mills theory through a Langevin equation 
[4] [5] . For reviews of stochastic quantization, see [6] [7] . The 5-dimensional theory will be 
a supersymmetric theory of the topological type, that we can express as a path integral 
over 5-dimensional gauge fields for the continuum perturbative formulation, or as a lattice 
gauge theory which now depends on fields of a topological field theory. In one of the lattice 
formulations that we will present, the stochastic time is discretized. In the 5-dimensional 
critical limit controlled by — 0, both the physical Euclidean correlation lengths and the 
dynamical Monte Carlo auto-correlation times diverge simultaneously. This allows new 
algorithms. 
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The link between topological field theory and stochastic quantization was first ob- 
served in [8] . The general idea is that the D + 1-dimensional supersymmetric formulation 
of a D-dimensional quantum field theory that is obtained by expressing the field quanti- 
zation by a Langevin equation, involves the fields of a topological field theory in D + 1 
dimensions. Then, supersymmetric cancellations wash away the detail of the theory in the 
bulk, while the relevant aspects of the physical theory are retained in some boundary of 
the D + 1-dimensional space. 

Moreover an extra gauge field component is needed to enforce the 5-dimensional gauge 
symmetry in the stochastic framework [9]. The latter can be identified as a potential 
for a drift force along the gauge orbits of the original D-dimensional gauge theory, as 
independently observed in [10] [9] [11]. A gauge-fixing drift force for stochastic quantization 
was actually introduced originally in [12] , as a function of the 4-dimensional gauge fields yl^. 
But the point is that, by promoting the potential for this force to a 5th field component, 
and by functionally integrating over all its possible values, [9], not only does one not alter 
the physical physical sector of the theory, but one softens the gauge condition of the D- 
dimensional quantum field theory. This actually gives the desired result that the Gribov 
ambiguity losses is relevance as an obstacle for the gauge fixing of the theory. 

By postulating that QCD4 should be considered from a 5-dimensional point of view, 
we will emphasize its relationship to a topological field theory which establishes a pleasant 
geometrical framework, both in the continuum and lattice descriptions. 

Renormalizability of stochastically quantized scalar 0^ theory was demonstrated in 
[13] using the BRST operator which encodes the supersymmetry of the stochastic process. 
Renormalizability of the 5-dimensional formulation of gauge theory was demonstrated in 
[14] using the BRST operator for gauge invariance supplemented by graphical arguments. 
In the present article we demonstrate renormalizability using the complete BRST operator 
that encodes both the supersymmetry of the stochastic process and the gauge symmetry 
of the theory, consistent with parabolic propagation of all ghost fields. 

The disappearance of the Gribov ambiguity comes as an immediate consequence of 
being in five dimensions, in the context of renormalizable gauges (adapted to the power 
counting of this dimension). We will check this result by the verification that the ghost 
propagators cannot have zero modes. Alternatively, the mathematically oriented reader 
will notice that the argument of Singer [2] for having a Gribov copy phenomenon just 
disapears because the stochastic processes is defined on an infinite interval, even if we have 
a compact 4-dimensional physical space. 
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The consistency of our approach is ensured by the topological invariance, defined 
modulo ordinary gauge symmetry. We will also indicate that instantons in four dimensions 
are replaced by 5-dimensional solitons, but we leave open the question of anomalies, which 
is presumably connected to interesting 5-dimensional topological questions. 

To make contact to non-perturbative physics, we will formulate the fields and the 
symmetry in a 5-dimensional lattice formulation. Not only this offers a new point of view 
for the topological symmetries, but this allows us to obtain a discrete lattice formulation 
with a BRST invariant gauge- fixing, offering thereby a new understanding of the continuous 
limit. It also gives a concrete definition of non-perturbative physics. 

In sees. 2-5 we present the continuum 5-dimensional topological formulation. In sec. 
6 we describe the lattice stochastic formulation including the lattice Fokker-Planck and 
Langevin equations, in which the 5th time t — is identified with the number of Monte 
Carlo sweeps. We describe a new Monte Carlo algorithm based on the 5-dimensional 
critical limit which controls both physical correlation lengths and dynamical auto- 
correlation times, with time-step in lattice units given by e ~ ^^^^^ in gluodynamics. 
We conjecture that the physical spectrum may be given by the eigenvalues of the Fokker- 
Planck hamiltonian. In sec. 7 we present a 5-dimensional topological lattice formulation. 
In sec. 8 we indicate how our approach extends to fermions, and we conclude with some 
speculation on future developments. Readers who are primarily interested in Monte Carlo 
calculations may read sec. 6 independently. 

2. The 5-dimensional continuum action 

2.1. The Langevin equation 

The Langevin equation proposed by Parisi and Wu was [4] : 

dr,A^ = DxFx^ + V,., (2.1) 

where ry^ is a Gaussian white noise. In the stochastic approach, the correlation functions 
are functionals of An{x'^,x^) that are the solutions of this equation. The 4-dimensional 
Green functions are obtained in the limit where all arguments in the stochastic time are 
equal, and x^ — > oo. Since we consider a stochastic process, we can as well take initial 
conditions at x^ ~ — oo, and define the Green functions at any given fixed time, for instance 
x^ = 0. The convergence toward an equilibrium distribution relies on ergodicity theorems 



3 



that accord with physical intuition for which the relaxation of a gas to the state of maximal 
entropy is obvious, whatever of the initial distribution of its constituents. 

A diflBculty with this Langevin equation is that it provides no restoring force along 
the gauge orbits, because the Euclidean action S{A) is gauge- invariant. Consequently 
the probability escapes to infinity along the gauge orbits, and there is no normalizable 
equilibrium probability distribution. This may be remedied by modifying (2.1) by the 
introduction of a gauge-fixing term that is tangent to the gauge orbit [12] , 

d5A^ = DxFx^ + D^v[A]+r]^. (2.2) 

The gauge- fixing term [A] has the form of an infinitesimal gauge transformation, and 
consequently it has no effect on the expectation-value of any gauge-invariant observable. 
The infinitesimal generator t;[A](x,t) may, in principle, be completely arbitrary, and the 
issue of a "correct" gauge fixing, without Gribov copies does not arise. 

The last equation may be written in a way which respects gauge invariance in five 
dimensions. The idea is to identify v[A] with an independent 5th field component 

v = A5{x,t) (2.3) 

which eventually be fixed in the functional integral by a gauge-fixing term in the action. 
Then the gauge-fixed Langevin equation (2.2) may be written as: 

- DxFx^ = r/^ (2.4) 

where = d^A/^ — D^j^A^. The relation 

A, = a-^dxAx . (2.5) 

will be imposed in the context of topological field theory by a 5-dimensional gauge condi- 
tion. 

We would like to obtain a functional integral representation for the Langevin equa- 
tion (2.4). This is a non-trivial problem because there is a gauge invariance, and the 
determinant of the map which connects the ^'s and the noise ry's has longitudinal zero 
modes. Instead of attempting to solve these questions by manipulations on the Langevin 
equation, we will pass directly to a 5-dimensional topological quantum field theory. Later 
with lattice regularization we shall verify the consistency of the 5-dimensional formulation, 
including the absence of Gribov copies, and we shall shall show that after elimination of 
auxiliary fields present in this theory, one recovers the Langevin equation and the equiva- 
lent Fokker-Planck equation. 



4 



2.2. The 5- dimensional BRST symmetry 



The expression of the complete symmetry of the theory boils down to the knowledge of 
a BRST operator which encodes at once the supersymmetry of the stochastic process and 
the gauge symmetry of the theory. This operator is that of the topological symmetry of a 
gauge field in five dimensions, defined modulo gauge transformations. There is not much 
choice for exhibiting such a BRST operator, and the solution was given in [9]. It involves 
a fifth dimensional component for the gauge field, that we denote as A^. Unavoidably, 
the vector fermion that is needed to express the Jacobian of the constraint (2.1) must be 
enlarged into a five dimensional vector (^P^, ^'5). If we denote by c the ordinary Faddeev- 
Popov ghost and $ its ghost of ghost, the combined BRST operator for the stochastic 
supersymmetry and the gauge symmetry follows therefore from the following geometrical 
equation [9]: 

is + d){A + c) + -\A + c,A + c] = F + "^f^dx^ + "ifr^dx^ + $ 

2 ^ (2.6) 

{s + d){F + fj,dx^ + -^5dx^ + ^) = -[A + c, F + fj,dx^ + -^5dx^ + ^] 

The second equation is the Bianchi identity of the first, thus one has = on all fields 
which ensures the consistency of this symmetry which mixes the Yang-Mills symmetry 
with the supersymmetry of the stochastic process. Here c, '^^ and are fermi ghost 
fields with ghost number A^^ = 1 and $ is a bose ghost field with ghost number Ng = 2. 
This equation gives, after expansion in ghost number: 

sA^ = + D^c 

SA5 = *5 + D5C 

s*^ = - [c, 

(2-7) 

S*5 = -i^5*-[c,*5] 

sc = $ - ^[c,c] 

s$ = -[c, . 

It is identical to the topological BRST symmetry for a Yang-Mills field in five dimensions. 



We need anti-ghosts and lagrange multipliers, and their s-transformations. They are: 



s*^ = b^- [c, 
s*5 = h- [c, *5] 



sfes = -[c, 65] + [$,^'5] 

s§ = ^ - [c, l>] 



(2.8) 



sv=-[c,v] + [*,^] • 



All fields in (2.7) and (2.8) are valued in the same Lie algebra representation as A. In 
the lattice formulation of this symmetry, it will be advantageous to redefine the anti- 
ghost and lagrange multiplier fields to eliminate the dependance on c in (2.8), that is, 
= b', sb' = 0, = b', sb' = 0, = f}', and sfj' = 0. 

There are five degrees of freedom for the choice of a dynamics, one for each component 
of A, with Lagrange multipliers fields b^ and 65. One of them will serve to gauge-fix the 
ordinary 5-dimensional gauge invariance, with the condition (2.5); the other four will 
enforce the Langevin equation (2.4), following the now standard methods of TQFTs. The 
gauge fixing of the longitudinal modes in ^' will use 77 as a Lagrange multiplier. 

2.3. The action 

The complete action must be s-invariant. We will gauge fix the 5-dimensional invari- 
ant: 



This term is invariant under any local shift in Af^ and A5, since the Lagrangian can locally 



conditions for the variations, and eventually on the ^^'s. The choice of this term may appear 
as quite intuitive: it enforces the idea that the stochastic time is unobservable, since it is 
independent of the metric component (755, and it is compatible with Yang-Mills invariance. 
We will use an s— exact term for fixing all the invariances and introducing all relevant 
drift forces of stochastic quantization. Remarkably, the topological BRST symmetry just 
described is precisely what is needed to do so, that is, to represent a gauge theory in 
stochastic quantization by a functional integral with a high degree of symmetry [9]. 




(2.9) 



be written as TYd^{F^FI^)-\-\d^{F^Fl}). Of course this means that we use special boundary 
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To impose all relevant gauge conditions, we take the BRST-exact action: 

+ $(a'*5-i^A*A + f[r7,$]) (^-^^^ 
+ *^(a^5-5A^A)). 

The consistency of this choice of this gauge function will be further justified by showing that 
it gives a perturbatively renormalizable theory. Here a, a' and f3 are arbitrary parameters. 
When j3 is non-zero, quartic interactions are introduced. Of course /3 = gives simpler 
expressions, but nothing forbids /? 7^ 0, (/5 = might be a stable fixed point under 
renormalization). [Alternatively one may make a linear gauge choice by setting /? = 
and taking d\^\ instead of Dx^x] The first two terms are invariant under 5-dimensional 
gauge transformations. The first term concentrates the path integral around the solutions 
of (2.4), modulo ordinary gauge transformations, while the second term fixes in a gauge- 
covariant way the internal gauge invariance of that one detects in the BRST variation 
s'^ ^ = D^<^ + .... The third term fixes the gauge invariance for A^,A^. As we will see 
in section 3, all these terms are essentially determined by symmetry and power counting 
requirements. 

We may introduce a background Yang-Mills symmetry, which transform all fields, but 
A, in the adjoint representation, while A transform as a gauge field. The shortest way to 
represent such a symmetry, is to define it through a background BRST symmetry, with 
generator a and background ghost uj (which do not appear in the action). Since s and cr 
must anticommute, and = 0, we easily find that we obtain both the action of s and a 
by extending (2.6) into: 

is + a + d)iA + c + ijo) + ]-\A + c + u,A + c + ijo]=F + r,dx^ + "^^dx^ + $ 

2' J A* ^2j;L) 

(s + cr + d)(F + -^^dx^ + -^r^dx^ + $) = -[A + c + a;, F + ^dx^ + ^sda;^ + 

When one expands the latter equation one must assign a new ghost number to a and a;, 
which is independent of that carried by s and the propagating fields. This determines 
all transformations of s and a on all fields, with = = sa -\- as = 0. One then 
observes that the first two terms in (2.10) are not only s- invariant as s-exact terms, but 
they are also cr- invariant, without being cr-exact. These terms involve gauge-covariant 
gauge functions for A and ^ in five dimensions. The aim of the renormalization proof. 
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using power counting, will be to show that these terms, up to multiplicative rescalings, are 
the only ones which satisfy both s and cr-invariances. The third term is not cr-invariant, 
and it will be necessary to prove that when the parameter a changes, this does not affect 
the sector of the theory defined by s and a invariances. 

To investigate the properties of our action, we must expand (2.10), which is a simple 
exercise. One gets: 

1 = J dx^dx^ TV( - h^b^ + b^{F^^ - DxFx^) 

+ (65- [C, ^'5])M5-9^AJ 

- *^(a(*5 + D^c) - d^^^ - d^D^c)). 

(2.12) 

In order to see the dynamical content of this action, we identify and eliminate the 
auxiliary fields. The equation of motion of 65 gives back the gauge-fixing condition that 
we wish to impose: 

aAs = dxAx. (2.13) 
Moreover the equation of motion of 6^ has the same form as the Langevin equation, 

F5^-DxFxf, = b^. (2.14) 
We also eliminate and 77 and we get, 

a'vE'5 = I^M*p-f[r?,$] (2-15) 

a'fi = a^5 + D^^^. (2.16) 

Thus A5 can be expressed as a function of A^, and \l/5 can be expressed as a fimction of 
This is the key of a gauge-fixing that does not suffer from the Gribov ambiguity, and 
the gauge condition aA^ = dxAx combines the virtues of axial and Laudau gauges. By 
integrating over all values of A^, the Faddeev-Popov zero modes that one encounters in the 
genuine 4-dimensional theory just disappear. It is tempting to compare the integration 
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over A5 to the integration over moduli that also solves the problem of zero modes of 
reparametrization ghosts in string theory. 

We will give much more details about the elimination of Gribov copies in sec. 7, 
although a Fourier transformation on the variable could help understanding directly 
how the multiple intersections of gauge orbit that occur in four dimensions could be split 
on different orbits in five dimensions. 

Finally, after integrating out 6^ and 65, we get: 

1 + J dx^'dx'' TriF^^.DxF^") - J dx^dx^ Tr( - ^F^^F''^ - ^^D^F^^DpFP^" 

+ $[*^, ^z^]) + + ^[$, $][$, $] ^2.17) 

One has = a~^d\A\ in this action. 

To derive Feynman rules, we examine the part of the resulting action that is quadratic 
in the fields. It is given by 

Io = Ia + U + U + Ic (2.18) 
lA=j dx^dx'TT{^{d,A^f + lidxdxA^f + l(a-2 - l){d^dxAxf) 
I^=J dx^dx^TT{^^{-d5^^ + dxdx^^ + (a'-i - l)d^dx^x)) 
U= J dx'^dx^TT{^-a'd5^ + d^d^^)) 
Ic= J dx''dx^W^{ - ad^c + d^d/^c) 

where we have dropped exact derivatives in Ia- 

An important observation is that all ghost-antighost pairs, and 
(0,^5) have a free action which is parabolic: it is first order in ^5, and the correspond- 
ing matrix of spatial derivatives is a negative operator. Consequently all ghost propa- 
gators are retarded, D{x,t) = for t = 0^5 < 0. For example for a = 1, all ghost- 
antighost pairs have the momentum-space propagator (iuj + k^)~^, with fourier transform 
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(2.19) 



D{x,t) = ^(t)(47rt)~^exp(— a;^/4t), where 9{t) is the step function. Moreover ghost num- 
ber is conserved, and an arrow which represents the flow of ghost charge may be assigned 
to each ghost hne: it points from the past to the future. On the other hand, the action of 
the A- field is second order in 85,30 A fines freely move forward and backward in time. (Its 
momentum-space propagator is ^c<;^ + (/c^)^ j = (2k'^)~^ (j^iu + k'^)~^ + {—iu + k'^)~^^ .) 
Ghost number is also conserved at every vertex. Consequently in a diagrammatic expan- 
sion, every ghost line may be followed from the point where it enters a diagram from the 
past, as it moves monotonically from the past to the future, until it exits into the future. 
(The transformation of a $ into two ^^'s and back is allowed.) It follows that in a diagram 
containing only external ^-lines, there can be no closed ghost loops apart from tadples. 
With dimensional regularization the tadpole diagrams vanish, and in this representation of 
gauge theory, ghost diagrams provide accounting checks on the renormalization constants 
that are expressed in Ward identities, as usual, but they do not appear in the expansion of 
correlation functions the ^-field. With lattice regularization, the tadpole diagrams yield 
the famous Ito term. 

Another way to see that the ghost determinants have no Gribov problems is to note 
that the initial-value problem for a parabolic equation has a unique solution. In the 
section devoted to lattice regularization, we will indeed integrate out the ghosts exactly, 
and explicitly verify the absence of Gribov copies. 

We now turn to the question of understanding the renormalization of this action. 

3. Perturbative renormalization 

The action (2.12) is renormalizable by power counting. Indeed, the form of the propa- 
gators implies that the canonical dimension of all fields with an index 5 is equal to 2, while 
other fields have dimensions 1. It follows that we start from an action where the coupling 
constant has dimension 0, and we can apply the general result that it can be renormalized 
order by order in perturbation theory by a finite number of counter-terms. 

The main question is of course to investigate the structure of the counter-terms. We 
will only sketch the demonstration that this action is actually stable under renormaliza- 
tion, which means that the generating functional of Green functions can be defined, order 
by order in perturbation theory, while satisfying the same Ward identities as the action 
(2.10), corresponding to its invariance under s-symmetry, S'0(4)-symmetry, ghost number 
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conservation, and to the choice of a hnear condition aA^ — df^An in the sector which violates 
the cr-in variance. 

The action (2.12) that determines perturbation theory can indeed be be spht in two 
parts: 

/ ~/(^p, ^5, *5, ^, V) 

( r \ (3.1) 

The first term, where the fields A^ and must be understood as independent fields, 
(before eliminating 65), is s- and cr- invariant, while the last part is only s- invariant. In 
order to achieve the algebraic set-up which ensures perturbative renormalizability while 
maintaining the form of the action (up to so-called multiplicative renormalization) , it is 
enough to remark that one can write the second term as: 

=h{aA5 - d^A^) - *^(aL»5 - d^'D^y ^^'^^ 

where c' = c+{ad^ — D^d^)~^{a^^—d^^^). The change of variable c ^ c' is pcrturbativcly 
well defined, with a trivial Jacobian, because the operator ad^ — D^dfj, is parabolic. (This 
will be explained in detail in section 7, devoted to the lattice regularization.) 

If we now introduce an ordinary BRST symmetry s', with c' as a Faddeev-Popov 
ghost, we are almost in the ordinary situation of the renormalization, except that we 
have more fields. We have that s' acts on A^^A^, ^f^, 6^, ^^5, ^, fj as ordinary gauge 
transformations with the parameter equal to c', s'c' = —c'c', = 65 and = 0. The 
action (3.2) can be written as: 

/ ^i{A^, A5, t^, bf,, *5, ^, v) 

(3.3) 

+ s'{^\aA5-d^A^)) 

The first term is s'-invariant as a consequence of its cr-invariance: it corresponds to the 
cohomology with ghost number zero of s' , restricted to the terms with the correct power 
counting (6 with our power assignmant), while the second term is s'-exact, and will remain 
stable order by order in perturbation theory because of the linear gauge condition aA^ — 
d^A^. Notice that / depends on a', but ^^-dependent observables are independent of this 
parameter, as will be shown in sec. 7. 

The complete proof that the theory can be renormalized multiplicatively is just a slight 
amplification of the standard proof of the renormalizability of the 4-dimensional theory in 
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a linear gauge. It involves introducting sources for all s' variations, as well as using the 
equation of motion of the antighost as a Ward identity. When the Ward identity of the 
s'-invariance is combined with that of the topological s-invariance, using locality, one can 
prove from purely algebraic considerations that the counter-terms can be built order by 
order in perturbation theory such that, when they are added to the starting action (2.10), 
they provide an identical action, up to mere multiplicative renormalization constants for 
the fields and parameters. Essential is in this proof is the use of a linear gauge function 
aA^ - d^A^. 

A notable feature of the renormalization is that besides the multiplicative renormal- 
ization of the fields and coupling constants there is also a renormalization of the stochastic 
time 

t = Zttr . (3.4) 

4. Observables and gauge invariance 

It is now clear how observables should be defined: they are the cohomology with ghost 
number zero of s'; moreover, they must be computed at a same value of x^. Using the 
translation invariance under x^, and provided that the stochastic process starts at x^ = 
—00, the physical theory can be defined in any given slice of the 5-dimensional manifold. 
This points out the relevance of the boundary term (2.9), which is independent of the metric 
component g^^. Perturbatively, the equal-time gauge- invariant correlation functions are 
guaranteed to agree with those calculated in 4-dimensional theory by the usual Faddeev- 
Popov method. On the other hand, gauge non-invariant correlation functions are not 
expected to agree within the framework of local renormalizable 5-dimensional gauges that 
we consider in this paper ^. 

There are also Green functions that can be computed with fields at different values 
of x^, as well as the Green functions that involves the topological ghosts. We leave their 
interpretation as an interesting open question, although they are not directly relevant to 
physics. 

Once the renormalization has been done properly, one can integrate out all the ghost 
fields of the type \E', $ and c, provided that one considers observables which depend on 

^ To recover the local 4-dimensional Faddeev-Popov distribution, one needs a non-local 5- 
dimensional gauge-fixing on A5 [15]. 
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A/^ only. This follows from the parabolic or retarded behavior of all ghost propagators 
in the 5-dimensional theory. It will be demonstrated by explicit calculation with lattice 
regularization in section 7 that, when the ghost fields are integrated out, only the tadpole 
diagram survives so the ghost determinent is trivial and contributes only a local term to 
the effective action which is in fact the famous Ito term. This result, which explains the 
disappearance of the problem of Gribov copies, has an immediate application: mean values 
of observables are independent of the parameter a'. It should be clear however that the 
ghosts and ghosts of ghosts are nevertheless necessary to unveil the topological and gauge 
properties of the theory and to control its Ward identities. 

The other question is obviously that of the independence of physical expectation values 
of the gauge parameter a which appears in the gauge function aAr, — d^A^. To prove this, 
one observes that the observables are s'-invariant without being s'-exact, while the a- 
dependence is through an s'-exact term, so the standard BRST method based on Ward 
identities applies. One can also prove directly the non-renormalization of the s'-exact term 



We mention the possibility of introducing the interpolating gauges used in [16] for 
defining the "physical" Coulomb gauge as the limit of a renormalizable gauge. 

Finally, one must of course verify that no anomaly in the s- and cr-symmetries can 
occur. The absence of s-anomalies is quite obvious if one examines the consistency condi- 
tions for a topological symmetry; for the background Yang-Mills symmetry, no anomaly 
is expected, since we are in the case of a pure Yang-Mills theory, without 4-dimensional 
chiral fermion. The inclusion of spinors is a most interesting question, on which we will 
comment in the last section. 

5. Non-perturbative aspects 

In the previous section, we have verified that perturbation theory is as well defined in 
the 5-dimensional approach as it is in the 4-dimensional one. 

5.1. Global properties of gauge-fixing 

The gauge choice made here has the global property of a restoring force because it is 
derivable from a "minimizing" functional. Consider the functional 



in (3.3). 




(5.1) 
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which is proportional to the Hilbert square norm. Here a > is a gauge parameter. (More 
generally, one can take J-'[A] = J d^xAxax^Afj,, where a is a strictly positive symmetric 
matrix that characterizes a class of interpolating gauges [16].) The gauge condition used 
here, aA^ = d^A^, may be expressed as 

A^ix) = -G{x)T, (5.2) 

where 

is the generator of local gauge transformations. It satisfies the Lie algebra commutation 
relations 

[C^ix), G\y)] = 5{x - y)r''^G%x) . (5.4) 
Indeed, one can verify that 

Gix)J^=-a-^d^A^{x). (5.5) 

To see that the gauge-fixing force is globally restoring, consider the fiow defined by 
the gauge-fixing force alone 

d^A^ = D^v = -D^GJ^, (5.6) 

in the Langcvin equation (2.2). Under this fiow the minimizing functional J-' decreases 
monotonically, since 



M (5.7) 



= -{GT, GT) < 0. 



A consequence of our gauge choice is that the Langevin equation (2.4) is parabolic. 
Indeed, let A^ = a~^dfj,A^ be substituted into this equation, and consider the highest 
derivatives on the right-hand side, 

riA 

= dxdxA^ + (a-i - l)d^dxAx + . . . , (5.8) 

where the dots represent lower order derivatives. One sees that for a > 0, the operator 
that appears on the right is negative, which assures convergence at large t. 
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5.2. How the Gribov problem is solved 

Let us now discuss the Gribov question in more detail, and explain how our use of the 
gauge function aA^ — d^A^ escapes the ambiguity. 

The origin of this ambiguity in the 4-dimensional formulation, as given by Singer [2], 
was the compactness of the (Euclidian) space. It does not apply anymore, since the 5- 
dimensional formulation is necessarily non-compact: the system requires an infinite amount 
of time to relax to equilibrium. In our formulation the stochastic time runs from —oo to 
oo, and one evaluates the Green function at an arbitrary intermediate time, say = 0. 

Then, there is the explicit argument, which we already mentioned, of considering all 
ghost propagators, whether they are the and $ topological ghosts or the Faddeev-Popov 
ghost: by taking the suitable boundary conditions at = — oo, no zero modes can occur 
since we have parabolic propagation. It is a convincing argument, and we will check it in 
great detail in section 7, in the framework of the lattice regularization of our 5-dimensional 
formulation. We will prove that the integration on all ghosts does not lead to zero modes 
that would make ambiguous their elimination. In our opinion, this a concrete verification 
of the assertion that the Gribov problem is solved. 

To see intuitively where the gauge-fixing force concentrates the probability, observe 
from (5.7) that the gauge fixing force is in equilibrium only where G{x)!F = 0. This 
equilibrium may be stable or unstable, but the gauge fixing force drives the system toward 
stable equilibrium only, namely where the second variation of J-'[A] under infinitesimal 
gauge transformation is a positive matrix. This is equivalent to the condition that the 
Faddeev-Popov matrix 

M'^\x,y) = G%x)G\y)J' (5.9) 

has positive eigenvalues. With GJ^ = — a~^5^^^, the matrix M and the condition for 
stable equilibrium are given by 

M = -a-^d^D^ > 0. (5.10) 

Thus the conditions for stable equilibrium are (a) transversality of the vector potential, 
d^A/j, = and (b) positivity of the Faddeev-Popov operator, —dfj^D^ > 0. These two 
conditions define the Gribov region, a region which Gribov suggested (not quite correctly) 
was free of copies. As long as the gauge-parameter a > is finite, the gauge-fixing is 
soft in the sense that the weight on each gauge orbit has a spread that is centered on the 
Gribov region. As a approaches 0, the equilibrium probability gets concentrated close to 
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the Gribov region. The gauge-fixing force drives the system by steepest descent along a 
gauge orbit toward a relative or absolute minimum of the functional Tj^j [g] = J^l^A] , so 
Gribov copies that are unstable equilibria (saddle points of I[a] [g] ) are avoided in the limit 
a — > 0. However the relative minima are in fact Gribov copies of the absolute minimum, so 
in the limit a — > this gauge-fixing distributes the probability in some way among these 
relative minimum, possibly even entirely at the absolute minimum. However the validity 
of the present approach in no way depends on taking the limit a ^ 0. 

It is interesting to ask how instantons play a role in the 5-dimensional presentation, for 
instance in a semi-classical approximation. Their interpretation is actually quite simple. 
The bosonic part of our action is a sum of squares, namely 1-^5^ P + \DxFx^\'^. It obvi- 
ously gives and an absolute minimum when it vanishes. This occurs for x^-independent 
solutions, with A5 = d^^A^ = and satisfying D\F\^ = 0. This means that the classical 
4-dimensional instanton solutions also minimize the 5-dimensional path integral: they are 
the solitons of the 5-dimensional theory. 

We are therefore in a situation where, from various points of view, the path integral 
seems well defined. We can now seriously consider its evaluation for "large" field config- 
urations, since the question of duplicating orbits does not occur. It is however clear that 
the space of connections is not appropriate for defining a meaningful measure, and we will 
shortly look at the lattice formulation of this extended version of the Yang-Mill theory, 
keeping in mind that the use of the fifth time allows a more complete presentation of the 
theory. 



6. Stochastic lattice gauge theory and new Monte Carlo algorithm 

6.1. Lattice gauge theory with discrete time 

Numerical calculations in lattice gauge theory are efi^ected by Monte Carlo methods 
which rely on the simulation of a stochastic process. The stochastic process in question is 
specified by a matrix T{Ut-\-ii Ut) > of transition probabilities Ut — > ?7t+i, on configura- 
tions of "horizontal" link variables, Ut — {Ux,t,ix} and Ut+i — {Ux,t+i,y}, where = 1, ...4. 
Consider a transition probability of the form 

r(t/,+i,C/0 =^exp{-/3 J]Relv(/- ) } , (g.l) 

where /? is a positive parameter. Here V^,t,ju = V^,ju(?7t) is a group element that depends 
on the horizontal configuration at time t that will be specified shortly. For appropriately 
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chosen normalization constant this expression satisfies the requirement that the sum 
of transition probabihties out of any state is unity 

/ n dU,,t+^,^T{Ut+u Ut) = l. (6.2) 

This follows from invariance of the Haar measure under translation on the group, 

dUa;,t+i,i^ = dU'a;,t+i,i^, where t/'x,t+i,M = yx,t,f^Ux,t+i,,^- 
We set 

Vx,t,^^ = U~l^ e-Kp{ef:,,t,n)U^,t,nU^+(i,t,5, (6.3) 

where e is another positive parameter, and fx,t,n is a lattice analog of the continuum drift 
force {DxFxfj)^°'^^ that is specified below. This gives 

x,ix 

(6.4) 

where the product of 4 Vs is the transporter around the plaquette in the (/i-5) plane, 
starting and ending at the site (a;,t). This expression is manifestly invariant under 5- 
dimensional local gauge transformations provided that fx,t,ix transforms like a site variable 
at X in the adjoint representation, 

fx,t,ix ^ ^ fx,t,ij, — {9x,t) fx,t,fj, 9x,t 1 (6.5) 

which will be the case. 

In the last expression for T, the variables C/a;,t,5 associated to time-like links serve 
merely to eff'ect a gauge transformation g^^t = Ux,t,5 on the variables Ux,t+i,fj, associated to 
the space-like links. Consequently the Ux,t,5 may be assigned arbitrarily without affecting 
expectation values of the observables which are the gauge-invariant functions of the space- 
like link variables at a fixed time. We shall shortly gauge-fix Ux,t,5 by a lattice analog of 
continuum gauge fixing implemented above, where by Ux,t,b will be expressed in terms of 
the Ux,t,n so that eq. (6.2) is satisfied. 

Recall that the 5-dimensional continuum theory contains no dimensionful constant, 
apart from the cut-off, when engineering dimensions are assigned according to [t] = [x'^]. 
Consequently the discretized action should depend on the lattice spacings only through 
the ratio 



at 

e = 

a 



2 ' (6.6) 



s 
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which defines the parameter e. Here ag and at are the lattice spacings in the space and 
time directions respectively, with at being the "time" for a single sweep of the lattice in a 
Monte Carlo updating. We write 

U^,t,^ ~ exp(^oa.^^f^) 

(6.7) 

Ua,,t,5 ~ exp{goatAl%) , 

where go is the unrenormalized coupling constant, and "cont" designates continuum per- 
turbative variables. 

To verify this point, we estimate the quantities that appear in eqs. (6.1) and (6.3). 
The lattice drift force, specified below, is of order 

Ut,^ goaliD^F^^r^' . (6.8) 

This gives 



and we obtain 



(6.9) 



(6.10) 



exp(e/^,t,^) ~ exp (goagat {DxEx^y""^^^ 
Ux,t,f^U^+f.,t,5U-l+^^^U-l^ ~ exp ( - goasat (^5,^)'°"') , 

and finally 

ReTv(l -exp(e/a;,t,^) U^,t,^.U^+^l,t,5U-l+-L^^U-l^^ 

where a sum on fx and a is understood. We have 

YlT{Ut+i,Ut)^exp{-/3I) , (611) 
t 

where the five-dimensional discretized action is given by 

/ = ^ Re Tv(l - ex^p{efa,,t,^,)U^,t,^^U^+fi,t,5U~l+l^^U-l^^ . (^g_i2) 

X,t,IJ. 

We define 

/9 = ^ = (6.13) 



- \^t< E ((^5M - DxF^X''')' ' (6.14) 



which is the correct volume element, a^aj, and the correct normalization of the 5- 
dimensional action. 
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There remains to specify the lattice drift force fx^n and the gauge-fixing of Ux,t,5- For 
this purpose we shall establish a correspondence between lattice and continuum quantities. 
Let S = S{U) be the gauge-invariant 4-dimensional Euclidean lattice action normalized to 

S{U) ^ g'.S"^^-' = ^J d^x{F;j^°-'f , (6.15) 

which depends only on the horizontal link variables U = {Ux,fj,} at a fixed time. It may 
be the Wilson plaquette action 



5w = 2^Re Tr(/- t/p) 



(6.16) 



for pure gluodynamics, where the sum extends over all plaquettes p, or the effective action 
that results from integrating out the quark degrees of freedom, 

SiU) = Sw- g^Tr Hl^.D^ + m) , (6.17) 

where ^^Df^ + m represents the lattice Dirac operator of choice. Let the lattice color- 
"electric" field operator Ex,n,a be defined by 

Ex,ij.,a = {taUx,ij)al3-Q7Tj x , (6.18) 

where [taih] = f°'^^tc, and tT{tatb) = —{l/2)5ab- It has 4 components, A* = 1, ...4, because 
our perspective is 5-dimensional. With {Ux,^)ai3 ~ ^ap + {ta)ai3 9oci'sAx^^°^^ , we have the 
correspondence between lattice and continuum quantities, 

r ccont 

Ex,^,aS{U) ^ goal -— ^ = -goaUDxP^.r^' , (6.19) 
and we take for the lattice drift force 

fx,^^,a = -Ex,^,,aS{U) . (6.20) 

The lattice color-electric field operator defined by the remarkably simple expression 
(6.18) represents the color-electric flux carried by the link {x, ji). It satisfies the Lie algebra 
commutation relations 

\Ex^fj,^ai -E'j/,1/,6] — ^Xjy^HjV f Ex^fx^c i (6.21) 
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as well as 

^ J (6.22) 

[■Ex,ix,a-,Uy^i,'\ = — (^a;,j/<^/u,i/ Uyj^ta ■ 

We may visualize the color-electric field operator Ex,^.,a acting on S{U) as inserting ta into 
the plaquettes that contain the link (x, /u) on the left side of the matrix Ux,^- 
The color-electric field operator that operates by right multiplication of ta-, 

d 

E'x,ix,a = {Ux,fita)a/3^77-r x (6.23) 

satisfies the Lie algebra commutation relations with opposite sign 

[E x,^l,a■l E y,T/,b\ — ^x,y^fj.,v f E x,ix,c • (6.24) 

Since left and right multiplication commute, we have 

[Ex,^,,a,E'y,,,b]=0 ■ (6.25) 

With 

Ux,^ta = iUx,,,taU-l) Ux,^ = ObaiUx,^) hUx,^ , (6.26) 

where the real orthogonal matrices Oi,a{Ux,n) = Oab{U~^^) form the adjoint representation 
of the group, the two are related by 

E'x,n,a = Oi)a{Ux,^j)Ex,^j,,h ; (6.27) 

and satisfy 

^ x,ix,a 2-^^x,\x,a • (6.28) 
a a 

The generator of infinitesimal gauge transformations Gx,a is easily expressed in terms 
of these operators. A generic infinitesimal gauge transformation Ux,^ g~^Ux,^gx+jii 
with = 1 -I- a;^;, is given by 

^Ux,n = —<jJxUx,n + Ux nOJx+a 

(6.29) 

= {DiJ,'^)xUx,m ) 

where 

{D^Uj)x = Ux,^UJx+fi.U~^^ -^x= ta{Dfj,Uj)x,a ■ (6.30) 
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This defines the lattice gauge-covariant difference that corresponds to the continuum 
gauge-covariant derivative. It has a simple geometric meaning, being the difference between 
uj dit X + p, transported back to x, and a; at a;. This infinitesimal gauge transformation 
induces the following change in a generic function J-{U): 



(6.31) 



X 

This defines the generator of local gauge transformations Gx,a- From 

{Dfj,(jj)x,a = Oab{Ux,n) i^x+fi,b " ^a;,o (6.32) 

we obtain 

which is the total fiux of color-electric field leaving the site x. As in the continuum theory, 
it is the left-hand side of Gauss's law. One easily verifies the commutation relations 

[Gl.E'lj=5x-ix,yr"'E'l, (6.34) 

f/oo /^6i r fabc/^c 

[^x^^yi — '^x,y J ^y ■ 

Under local gauge transformation, E^ ^ transforms like a site variable at x in the adjoint 
representation, as required for the gauge-invariance of the transition matrix T{Ut+i,Ut). 
Moreover E'^ ^ transforms like a site variable in the adjoint representation at x + jj. 

We shall implement a lattice gauge fixing that corresponds to the continuum gauge 
fixing described above. We introduce a minimizing function, T{Ut), and we impose the 
gauge condition 

Ux,t,b = exp(^a;,t,5) = exp - eGx,t^{Ut)^ . (6.35) 

With this definition, eq. (6.2) is satisfied, because both Ux,t,5 and the gauge-fixing force 
fx,t,ij. are expressed in terms of variables that live only on the hyperplane t. For the 
minimizing function we take a lattice analog of J d'^x{A'^°^^)'^, for example, 

.F(f/) = 2a-^5^ReTr(l-f/,,^) 

""'^ (6.36) 
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where a > is a gauge parameter (not to be confused with the lattice units and at). 
This gives the gauge condition 



= —ea 



(6.37) 



6.2. Five- dimensional critical limit and new algorithm 



The renormahzabihty of the 5-dimensional continuum theory demonstrated above 
controls the critical limit of the 5-dimensional lattice theory just described. In particular 
it fixes the dependence of the parameter e = e{go) on go, and thereby provides a new 
accelerated Langevin algorithm, as we shall see. 

The probability for transition from configuration Ut into the volume element 
Ilx,nd'Ux,t+i,i^, where dUx,t+i,fj, is Haar measure, is given by 



Y[du^,t+i,^.r(Ut+uUt) . 



x,n 



(6.38) 



From eq. (6.4), and invariance of Haar measure under translation on the group, it follows 
that the probability for transition from the configuration Ut to the configuration 

Ux,t+l,^l = U'l^ exp{efa;,t,^)Ua;,t,i^Ua;+ii,t,5Wa;,t,^ (6.39) 

depends only on the group element Wx,t,iJ.j 

n dU,,t+i,^T{Ut+i,Ut) =Afll{dW,,t,^, exp ( - pRe Tr(/ - W,,t,^)^ }. (q^q) 

X,IJ, X,IJ, 

Thus to generate a new configuration it is sufficient to generate the last probability dis- 
tribution independently for each group element Wx,t,n, and also independently of the 
configuration Ut, and then calculate Ux,t+i,^i from eq. (6.39). (If one is not interested in 
gauge-fixing, one may set Ux,t,5 = 1- However the gauge-fixing defined above smooths out 
configurations, which may help to accelerate the inversion of the Dirac operator.) One 
recognizes the familiar Langevin algorithm for the SU(N) group [17] and [18], with Wx,t,n 
being white noise on the group. The new element is that we shall determine e = e{go) from 
the perturbative renormalization group. 
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We noted in our discussion of the continuum theory that the time is renormahzed 
according to t = Zttr- The parameter e = at/a^ rescales in the same way as the time t, 



e = Zter , (6.41) 

where is independent of the ultraviolet cut-off A = a~^. Consequently it satisfies the 
renormalization-group equation 

.dine ^ dlnZt 
dA dA 

^ Ptigo) (6-42) 

= -bt,o9o + Oig^) , 

where l3t{go) is a /3-function for the renormalization of the Monte Carlo time. The coeffi- 
cient bffi has been calculated for a theory without quarks, [19] and [20], and is independent 
of the gauge parameter a, 

It is convenient to change dependent variable from A to ^o? using 

, dine ^, .dine . .x 



where 



is the usual /9-function, 



P{9o) ^ A^ (6.45) 



/3(^o) = -hogl - higl + 0{gl) . (6.46) 



It has been verified to one-loop level for the 5-dimensional theory without quarks that 
/?(^o) is the same as in the Faddeev-Popov theory, [19]and [20]. For the theory without 
quarks we have 



which gives 



so 



c^lne ^ Ptjgo) 

dgo Pigo) ^g4g^ 

bt,o -1 N 
= ~^9o +O{go), 

= {hoglt^'^^''^ exp (0{gl)) C , (6.49) 
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where C is a constant of integration, and 

^ = (6.50) 

for the theory without quarks. The higher loop corrections are negligible in the critical 
limit go 0. This allows us to take for the purposes of numerical simulation 

e=(Mo)''°/^''°^ (6-51) 

and 

P = {g'oe)-' = ((^0')'+''°/^''°^ . (6.52) 
In a theory without quarks, this gives 

e = {hogl)-'"^' C 



(6.53) 



We note that (3 diverges in the critical limit ^0 — 0, but more slowly than I/qq. 
Moreover the time-step in lattice units diverges, e = a^/s^ —> 00 because 60 and 6o,t have 
opposite sign in a gauge theory. This provides a highly accelerated algorithm. However 
one may doubt whether this algorithm converges because high wave-number components, 
A; ~ aj-^ will overshoot and thus may appear to fluctuate wildly. However they are not of 
physical interest in the critical limit. On the other hand the modes that are of physical 
interest k ~ Aqcd will not overshoot. Indeed the real problem is likely to be that they 
evolve too slowly as usual, in which case the divergence of the step-size e{go) — > 00 as 
(7o — >^ is a great asset. Obviously the algorithm must be studied in practice before 
conclusions can be drawn. In order to apply the algorithm proposed here to a theory with 
quarks, the coefficient bt^o must also be calculated with quarks. 

If one takes C 00, one obtains a difi'erent, or non-existent, critical theory from the 
continuum theory discussed above. If one takes C ^ 0, one obtains the continuous-time 
stochastic theory discussed in the next section with a Euclidean spatial cut-off'. 

The 5-dimensional critical limit is approached in the limit 5^0 ~^ 0. The value of qq 
determines the physical correlation length as measured in Euclidean lattice units in the 
same way that it does in the Faddeev-Popov theory, 

1 \ 



^QhDibog'o)-'^/^'''^^ exp ( - ') , (6.54) 
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where Aqcd is the standard QCD mass scale, although this has been verified to one- loop 
level in the theory without quarks [19] and [20]. This formula and the above expression 
for €{go) fix the "time", at, for a single Monte Carlo sweep in terms of Aqcd, 

at = eCg 

a, = A5?,„C(6„92)<'"'>"-='>')/<-S' exp ( - ^) . '''''' 

In the asymptotic scaling region, auto-correlation times of physical observables are finite 
and independent of go and C when measured in units of Aq^j-,. The criteria for choosing 
the value of C are similar to those for ^ro- The auto-correlation time r for a physical 
observable should be large compared to one sweep-time, r >> at, but small compared the 
total running time T, r « T, so that the statistical uncertainty, of order N~^^'^, is small, 
where N = T/t. 

If it could be established that Monte Carlo auto-correlation times r of physical ob- 
servables are related to physical Euclidean correlation lengths L, by a universal relation 
T = KL^, where K is the same for all physical observables, then physical information such 
mass ratios, could be extracted directly from Monte-Carlo auto- correlation times. In this 
case the mass spectrum is given by the eigenvalues of the Fokker-Planck hamiltonian, as 
explained in sec. 6.3. 

We briefly mention some attractive features of the algorithm proposed here, which 
are those of the Langevin algorithm with a time-step in lattice units that diverges in the 
critical limit ^ro oo- There is no accept/reject criterion. Every updating Wx,t,fj, is 
accepted, even when the (effective) action is non-local as occurs with dynamical quarks. 
Only the derivative of the Euclidean action S appears in 

fx,t,^. = -E,,^.S = -E,,^Sw + g^Tr(^{^xDx + m)-^ E^^^{^,D,)) , (6.56) 

so one does not have to calculate the determinent of the lattice Dirac operator, but only 
its inverse. Moreover it is sufficient to calculate this inverse only once per sweep, because 
a whole sweep of the Euclidean lattice, described by eqs. (6.39) and (6.40), corresponds 
to a single hyperplane of the 5-dimensional lattice with action (6.12). 

Of course other algorithms may be invented that are based on other discretizations in 
which the 5-dimensional critical limit is also achived in the limit ~^ 0. 
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6.3. Lattice Fokker-Planck and Langevin equations 

One would be tempted to formulate a BRST version of the lattice gauge theory de- 
scrobed in the previous section, which has discrete stochastic time t and group-valued 
link variables Ux,t,5 asociated to the vertical (time-like) links. However Neuberger [21] has 
shown that lattice BRST gauge-fixing fails if the gauge-fixing function is continuous on the 
group manifold. Moreover Testa [22] has provided a counter-example which shows that if 
one tries to get around Neuberger's theorem by choosing a gauge-fixing function that is 
not continuous on the group manifold, then inconsistencies in the BRST procedure appear 
in the form of BRST-exact quantities sX with non-zero expectation-value < sX >^ 0. 
Therefore in order to comply with Neuberger's theorem we shall make t = continuous 
while keeping a;^ discrete, so Ax,5 is in the Lie algebra whereas Ux,fj, is in the Lie group, 
where = 1, ...4. Gauge-fixing will be of the form Ax^5 — Vx{U), so that a quantity in 
the algebra is fixed instead of a quantity in the group. An alternative approach to BRST 
gauge-fixing in lattice gauge theory has recently been provided by Baulieu and Schaden 
[23]. 

Before developing the BRST formulation of this 5-dimensional lattice theory, we 
first consider the stochastic process described by the lattice Fokker-Planck and Langevin 
equations with random variable which is a 4-dimensional Euclidean lattice configuration 
U = {Ux^fj,}, where Ux,^ G SU{N). Note that a Euclidean lattice configuration U is a 
point on a finite product of Lie-group manifolds and that a continuous stochastic process 
on Lie-group manifolds is the probabilistic analog of quantum mechanics on Lie-group 
manifolds. Both allow a well-defined path- integral formulation. 

The continuous-time stochastic process is described by the Fokker-Planck equation 

BP 

^ = -i^ppP , (6.57) 

where P = P{U,t) is a probability distribution, and the Fokker-Planck hamiltonian is 
given by 

t ^^'( a -2 a \ a -2 a a ^^'^^^ 

X 

where 

Ax,5 = -GxHU) , (6.59) 
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and D/j, is the lattice gauge-covariant difference, eq. (6.30). In Appendix A we show that 
this equation describes the hmit of the discrete-time stochastic process defined previous 
paragraph in which the time-step approaches zero, e = a^/a^ — > 0. 

In the previous paragraph it was conjectured that the auto-correlation times r of phys- 
ical observables might be related to Euclidean auto-correlation lengths L by the universal 
relation r = KL^, where K is the same for all observables. If this is true, then the mass 
spectrum is given by the eigenvalue equation 

. (6.60) 

But for the gauge-fixing "force", {D^A^)^^ in i^pp, the equilibrium solution, satisfying 
HppPeq = would be given by Peq = ^w{~9o^ cannot give the exact equilibrium 

solution for non-zero gauge-fixing force because it is not conservative. However with lattice 
regularization gauge-fixing may be dispensed with. Then a similarity transformation by 
exp[— (5f^^/2)S'] is sufficient to make the the Fokker-Planck hamiltonian hermitian 

^FP = - Et i^l^ - i9o'm{El,S)) (eI^ + {9o'mEl,s) ] , (6.61) 

and Hpp = H'pp. Moreover the physically relevant eigenfunctions are gauge- invariant, 
Gx^' = 0. The eigenvalue equation 

H^p^' = m^*', (6.62) 

is interesting. It is both Euclidean- and gauge-invariant. Moreover the ground-state or 
vacuum wave-function is known exactly, = exp[— (^^^/2)5'], so variational calculations 
of excited-state eigenvalues are possible. 

The corresponding Langevin equation is obtained from the e — > limit of the stochastic 
process defined in eqs. (6.39) and (6.40). It is given by 

Ki^U-j^ = go\D^A^), - g^^E,,^S + r/,,^, (6.63) 

where rjx,ij, is Gaussian white noise with two-point correlator 

< r]x,i^{^)r]y,iy{t) >= '^5x,y5^,,v5{s - t). (6.64) 
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The gauge choice Aj;^^ = —GxJ^{U) has the global property that for the flow deter- 
mined by the gauge-fixing force alone, 

d5Ux,^U-l = g^\D^A^%, (6.65) 

the minimizing function T decreases monotonically 

d^T < , (6.66) 

for we have 

3C ^ jAt CC CC 

For the minimizing function (6.36) , this drives all the Ux,iu, toward unity. 



(6.67) 



7. Five-dimensional topological lattice gauge theory 
7.1. Lattice BRST operator 

We determine the lattice BRST operator s which corresponds to the continuum BRST 
operator defined previously. It possesses the symmetries of the hypercubic lattice. In 
order to exhibit these symmetries, we shall, in this subsection only, use a symmetric 5- 
dimensional notation, whereby x = {.x^} for |U = 1, ...5 represents lattice sites (previously 
denoted (a;,t)), and we shall also denote link variables by Uxy € SU(N), with Uyx = 
{Uxy)~^, where (xy) is a pair of nearest neighbors. 

Consider a generic infinitesimal transformation Uxy, 

Uxy + SUxy = (1 + i^xy)Uxy = Uxy{l + io' xy) , (7.1) 

where uJxy and uj'xy = {Uxy)~^uJxyUxy are elements of the Lie algebra, Uxy = ta^^xy We 
write 6 = es, uJxy = ^^xy and u'xy = ^^'xy, where e is an infinitesimal Grassmann variable, 
and Qxy and Q'xySixe Lie algebra-valued Grassmann variables, with 

^ xy — (Uxy) ^xyUxy • ('^•^) 

This determines the action of s on Uxy, 

^Uxy — ^xyUxy — Uxy^ xy • ij^'^) 
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From the condition = 0, we obtain 

Silxy = illy 

(-)/ _ _o'^ 
SiL xy — xy • 

We will not write out explicitly further relations for fl'xy since they may be obtained for 
those of Qxy 

As in the continuum theory, we wish to distinguish the infinitesimal gauge transfor- 
mations from among all possible transformations of Uxy We pose Ux = ecx, where e is an 
infinitesimal Grassmann variable, and Cx is a Lie algebra-valued Grassmann variable, so 
the infinitesimal gauge transformation defined in the previous section reads 

^gUxy — ^( CxUxy ~l~ UxyCy^ ^ 

— ~l~ UxyCyUyx^Uxy 1 

where we have written Uyx = {Uxy)~^ ■ We separate this infinitesimal gauge transformation 
out of ^x,ni and write 

^xy — ^ xy ^x ~l~ UxyCyUyx , (7.6) 

which defines '^xy The action of s on Uxy now reads 

^Uxy — xy Cx)Uxy ~\~ UxyCy . (T.T) 

We identify the lattice site variable Cx with the corresponding continuum scalar ghost c{x), 
and similarly for ^xj ^ind assign to them the same transformation law as in the continuum 
case namely 

^Cx — ^x C„ 

(7.8) 

S^x = -Cx^x + ^xCx ■ 

The action of s on ^xy is obtained from 

xy — ^^xy ~l~ ^Cx ^(UxyCyUyx) • (7.9) 

This determines the action of s on all the lattice fields and ghosts. 

We write it in the notation in which it will be used, namely with lattice sites designated 
by {x,t), where x, is a 4- vector x = {x^} for /i = 1, ...4, and t = x^. "Horizontal" links, 
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corresponding to links in 4-dimensional Euclidean space-time, are designated by {x,t,fi), 
and "vertical" links by {x,t,5). In this notation s acts according to 

sUx,t,n = {^x,t,i^ + {D,j,c)x,t) Ua:,t,^i 
sU^,t,5 = (*x,t,5 + (-D^c)a;,t) Ux,t,5 

S^x,t,5 = -(i^5^)x,t - [Cx,t,^x,t,5] + i^x,t,5f 
SCx,t = ^x,t - {Cx,t) 

= -[Cx,t, ^x,t] , 

which satisfies = 0, and where the lattice gauge-covariant difference operator is 
defined in (6.30). Under gauge transformation '^x,t,i^ and '^x,t,5 transform like site variables 
located at (x, t). 

We shall shortly take t = X5 to be a continuous variable, while x remains the discrete 
label of a 4-dimensional lattice site. In this case, instead of the link variable Ux,t,5 in the 
Lie group, we require a variable Ax,t,5 in the Lie algebra. The BRST operator acts on 
Ax,t,5 and 'ifx,t,5 as in the continuum theory, 

sAx,t,5 = *a;,t,5 + D^Cx^t 

(7.11) 

S'^x,t,5 = -D5^x,t - [Cx,t,'^x,t,5\ , 

where is the continuum covariant derivative and s acts on the remaining variables as 
above. 

In terms of the lattice color-electric field operator Ex,^ and generator of local gauge 
transformations Gx obtained above, the explicit form of the lattice BRST operator acting 
on functions / = f{U, c, $) is given by sf = Qf, where Q is the differential operator 

Q = Y1 + CxGx + {-{D^.^)x - [Cx, + "^Uq^ 

"^'^ (7.12) 

+ i^x-ci)- [c^,$. 



dcx ' 'd^^ 

The action of s on the anti-ghosts and Lagrange multipliers is defined simply by 

S^x,t,ix = bx^t,n 
S^x,t,5 = K,t,5 

sbx,t,i^ = , , 

(7.13) 

^^x,t = fjx,t 

sVx,t = • 
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This is the same as in the continuum theory after the change of variable 

K,t,5 = Kt,5-[Cx,t,^x,t,5] (7.14) 
'n'x,t = ^x,t - [Cx,U ^x,t] ■ 

Note that ^x,u^x,uK,t,ti and ha;,t,5 are real variables. 

It is remarkable that the lattice BRST transformation requires only minimal change 
from the continuum BRST transformation, even though there is no local lattice curvature 
and no lattice Bianchi identity, whereas our starting point for the continuum transforma- 
tion was the curvature equation and the Bianchi identity (2.6). 

7.2. Five- dimensional topological lattice action 

The most powerful representation of the continuous stochastic process just defined is 
by path integral. We start with the 5-dimensional topological lattice action which we write 
by analogy with the 5-dimensional continuum action presented above, 

X 

+ ^xMx,^ + Gxr) ^^-^^^ 

where a sum on y and is understood. Here J-{U) is the minimizing function introduced 
above, and !F'{U) is the same or another minimizing function. 
We expand the above action and obtain 

/ = /i +/A5 +/*^ +/c + /* , (7.16) 

X 

I A, = ^o"' fdtJ2^{ ^"'-'5 (^-.5 + Gx-^)) (7.17) 

X 
X 

= -9o^ I dtY,T^{^x,^, { D^(^x,^ + {D^c)x) - £>^(*^,5 + {D^c)x) 

+ ( + {D.c)y )Ey,A,t.S{U) } ) 
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=g^^ f dtY,M^x{- {D5^)x - [Cx, *a;,5] 

(7.20) 

+ ( - {D^^)y - [Cy, *,,^] + ^l^) Ey,^G,J^' } ^ . 

The terms Ia^ and I^g serve to gauge fix A5 and ^^5. Indeed upon integration with respect 
to b'5 and 77' they impose the constraints, 

(7.21) 

7.3. Proof of gauge invariance and absence of Gribov copies. 

We must show that the expectation-values of physical observables, namely the Wil- 
son loops, are independent of the gauge-fixing functions JF and JF'. For this purpose we 
integrate out rj' and ^'5, which results in being assigned its gauge-fixed value. We next 
integrate out the bose ghosts $ and $, which appear only in 

Z^ = j d^d^exp{-U) . (7.22) 

The integral on $ yields 

Z^ = jl[d^t,.l[s(^-d^^a> + {mx + fx) , (7.23) 

where fx is independent of and the finite matrix L is defined by 

(L$)« ^ -[Ax,„ - G'yGlT' $^ , (7 24) 

y 

and we have used J2fj.i^t^)y,t-i = To evaluate we observe that the stochastic 
process starts at some initial time t = to, which may be taken to be to = "Oo. At the 
initial time, ^x{to), and the other fields ?7, c, are assigned some definite initial value 
that is not integrated over, for example $a;(^o) = 0. With this initial value, the differential 
equation 

^50, - iL(l>)x = fx (7.25) 
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is equivalent to the integral equation 

cl>Ut) = f^du[fl{u) + L'^liuWM)- (7.26) 

This equation possesses a unique solution (j)x{t). Thus unlike gauge-fixing on the group, or 
in a 4-dimensional covariant gauge, there are no Gribov copies. 

We shift $ by which cancels fx in the path integral for We also multiply as 
usual by the formally infinite normalization constant det^s. This gives 

Z^ = j \{d^x{t)\{5[^x{t) - {m.{t)) , (7.27) 

where £ = ^^"'^L is the integral operator 

{C^)l{t) = / du Cf^yit, u)^l{u), (7.28) 

with kernel 

Ctiit,u) = d{t-u)L-^^l{u) . (7.29) 

This gives 

= det"^(l - £) = exp (I - Tr ln(l - £)) = exp (TtC + , (7.30) 
Because jCx,y{t, u) is retarded, jC{t, tt) = for tt > t, only the tadpole term 

dt 

-J x,a 

dtJ^^Zit) 

-J x,a 



(7.31) 



survives in the last expansion. It is natural to assign the ambiguous expression ^(0) its 
mean value 1/2, which is consistent with other determinations (see Appendix B), and we 
obtain 

Z$=exp(-/ dtY^LZit)) , (7.32) 

•'^o x,a 



and so, by eq. (7.24), 

= exp ( - - / dtJ2 GIGIT') . (7.33) 



*o a;,o 
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We emphasize that we have made an exact evaluation of the ghost determinant, and 
obtained a purely local contribution to the effective action. 

In the same way we integrate out the Fermi ghosts and after setting to its 
gauge-fixed value. They only appear in the factor 

Z^^ = j d*^dt^exp(-7^J . (7.34) 

As before only the tadpole term survives and gives 

/I ^°° \ 
= exp (- / dtJ2iJ2 + GIGIJ^')) . (7.35) 

"^^o x,a n 

The dependence on JF' cancels in the product 

Z^^Z^ = exp(-/2) , (7.36) 

where 

^2^-2/ 5Z • (7.37) 

Thus physical observables are independent of the gauge-fixing function T' ^ as asserted. 

We next show that they are also independent of T. The integration on h^' is effected, 
which causes Ii to be changed to /i', given below. The variable of integration c may be 
translated, just as $ was translated above, to cancel the inhomogeneous term in 1^ which 
changes Ic to Ij given below. Putting these factors together we obtain 

[H (dU:,,^{t)dA^,5{t)dc^{t)d^xAt)dbx,5{t)) exp(-/') , (7.33) 



t,x 



where 



/' = // + /2 + /^^ + Ij 

h' = -l9o' [ dtJ^T^ {iD,U),,^U-l + E,,^SiU)y 



(7.39) 



and I2 and Ia^ are defined above. 

We next use standard BRST arguments to show that the expectation-value of any 
gauge-invariant observable is independent of the gauge-fixing function J^. Observe first 
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that Ii and I2 are gauge invariant. They are therefore invariant under a "httle" BRST 
operator s' that acts according to 

s'c, = -cl (7.40) 
s'^x,b = bx,5 s'bx,5 = , 

which satisfies s'^ = 0. Indeed, s' is the usual BRST operator with Cx ^x,5 and 
bx bx,5, and it acts on Ux,^ and Ax,5 hke an infinitesimal gauge transformation. Because 
Ii' and I2 are gauge- invariant they are also invariant under s', 

s'h' = 0; s'h = . (7.41) 

Moreover we may express the total action /' as 

/' = Ii' + /2 + s'J , (7.42) 

where 

J = g^^ fdtJ^T^ {^x,5 ( Ax,5 + GxT )) . (7.43) 

x 

Thus /' is s'-invariant, s'l' — 0, and moreover the gauge-fixing function JF appears only in 
the term s'J that is s'-exact. By standard arguments it follows that the expectation- value 
of any gauge-invariant observable is independent of J^, as asserted, which completes the 
proof. 

We next integrate out the remaining auxiliary variables 65,^5,^ and c to obtain a 
path- integral in terms of the actual random variables Ux,^{t) only. In eq. (7.38) , we 
integrate out bx,b{t) and Ax,5{t)^ which fixes Ax,^ = —GxT{U). The integral on the Fermi 
ghosts Cx{t) and ^x,b{t) namely 

= / ]][dcx{t)d^xAt)\ exp(-/c') , (7.44) 

t,x 

may be effected just like the integral on the bose ghosts ^x{t) and ^x{t) that was done 
explicitly above. Again only the tadpole term survives and gives 



Zc = exp(-/3) 

(7.45) 

'to 



1 r°° 

^3 = -2/ dtJ^G^G-J^ 

■^*o x,a 
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This yields the desired partition function of the stochastic process in terms of the random 
variables Ux,ix{t) only, 

Z = I Y[dU,,^{t) exp(-/") , (7.46) 

t,x 

where 

I" = 1^' + 1^ + 1^ (7.47) 

is a local action given by 

•'^o X (7.48) 

The last two terms are the famous Ito terms that are discussed in Appendix B, and we 
see that they are produced automatically by the 5-dimensional topological action when all 
auxiliary and ghost fields are integrated out. Thus the ghost determinant is well-defined, 
and moreover its explicit evaluation gives a local Ito term to the effective action which 
comes from the tadpole diagrams only. The Ito terms, which involve the second derivative 
of the action S and of the minimizing function T with an over-all minus sign, favor minima 
of S and T over maxima and saddle points. This is physically natural and transparent, so 
it is somewhat surprising that the tadpole diagrams, which give the Ito terms with lattice 
regularization, vanish with dimensional regularization [14] . 



8. Conclusion and Perspectives 

We have seen that the description of the 4-dimensional Yang-Mills theory in the 5- 
dimensional framework allows one to obtain a field theory description that does not suffers 
from the contradictions of the four dimensional description. From a physical point of 
view, we have introduced in a gauge-invariant way the additional time of the Monte Carlo 
description of the quantum theory and unified it with the Euclidean space-time coordinates. 
It is of course quite remarkable that this approach permits one to get a fully consistent 
lattice formulation of the theory with a natural limit toward a continuum formulation, 
while the Gribov ambiguity is just absent. Our work shows that this ambiguity is merely 
an artifact of a purely 4-dimensional description. It is only in five dimensions that one has 
a globally correct gauge fixing which is purely local. It thus appears that we can consider 
the 4-dimensional physical space as a slice in a 5-dimensional manifold, in which one can 
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express the Yang-Mills theory under the form of a particular topological field theory. A 
way to this this schematically is to observe that the 4-dimensional action can always be 
rewritten as the integral of a topological term in five dimensions: 



Thus, the definition of a quantum field theory amounts to the BRST invariant gauge-fixing 
of the action, a task that we have explained in detail above and which leads to the solutions 
of several problems of the purely 4-dimensional formalism. 

We briefly indicate how our method extends to fermions. To obtain the 5-dimensionsal 
fermionic action of a 4-dimensional Dirac spinor q, which automatically extends to a 5- 
dimensional spinor, it is now quite natural to introduce a topological BRST operator, and 
to enforce the relevant Langevin equation by a BRST-exact action. Call '^q and the 
commuting topological ghost and anti-ghost of the anti-commuting spinor g, and hq the 
anti-commuting Lagrange multiplier field. Then one has sq = '^q — cq, = — c\l/q + $5, 
and one may take s'^q = hq, shq = or s^q = bq — c'^q, shq = — c6g + $^g. In the Euclidean 
formulation the anti-commuting anti-quark spinor field q^ is an independent field, and it 
has corresponding commuting topological ghost and anti-ghost and and Lagrange 
multiplier field hq, with sq^ = — q^ c, = ^Jc — q^(^. For the anti-ghost one has 
again s*J = b\, sh\ = or = h\ + ^\c, sb\ = -h\c - *|$. AU these fields are 
functions of x and t. We take the topological action. 



which must be added to the 5-dimensional action that we have introduced in the pure Yang- 
Mills case. Here (7^!)^ + 'm)q = is the Dirac equation of motion of the 4-dimensional 
theory and X is a kernel [24]. One can formally prove, by mere Gaussian integrations, 
that the equilibrium distribution is independent of the choice of the kernel K [5]. Since we 
are in the context of a renormalizable theory, it is in fact necessary to introduce a kernel 
with canonical dimensional equal to 1, that is K = — 7'*!)^ -|- M. This is remarkable, since 
this coincides with the necessity of having a well-defined Langevin process for spinors, as 
discussed for instance in [5]. Quite interestingly, power counting would allows a term of 




(8.1) 




/ 



d^xdt {hliD^q - K[{YD^ + m)q + ahq] } 



(8.2) 
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the type F'^^F^j^ that is not proportional to the kernel: such term would contradict the 
stochastic interpretation of the 5-dimensional theory and could be the signal of an anomaly. 
For vector theories, such term should not be generated perturbatively. 

Beyond solving the ambiguities which occur in a purely 4-dimensional approach, our 
work might offer new perspectives. We have already mentioned in sec. 6 the existence of 
a new algorithm and the conjecture that the mass spectrum is given by the eigenvalues of 
the Fokker-Planck hamiltonian. 

One speculation is the question of understanding chiral fermions in a lattice formu- 
lation. The topological fermionic action (8.2) is interesting from this point of view. One 
may consider introducing a t-dependent mass term //(t) in (8.2). This might lead to an 
action of the type used by Kaplan and Neuberger to solve the chiral fermion problem for 
the lattice [25] [26] [27]. In their picture, one introduces a tower of additional fermions, 
which can be combined by Fourier transform into a 5-dimensional fermion. One notices 
that the equation of motion of (8.2) reproduces their 5-dimensional Dirac equation when 
fi{t) is a step function. Moreover, if one gauges the fermion conservation number in the 
fifth dimension, which gives an additional U{1) gauge field Ag_]3ei{x,t), then fj,{t) could be 
understood as a special value for the gauge field component ^sabeK^^O; similar to a kink 
that cannot be set globally to zero by local gauge-fixing. Therefore the choice of /x(t) 
could be determined by topological considerations. This issue of seeing whether one can 
understand the 4-dimensional space in which the stochastic process stops as a domain wall 
in which everything that remains from the fields of the topological multiplet of g is a chiral 
fermion seems to us to be great interest. 

The second perspective is that of a dual formulation which would now involve 2-form 
gauge fields. It is worth indicating the idea here. In the 5-dimensional framework, the 
dual of a one-form is a 2-form, since the sum of the degrees of the curvatures of A and its 
dual must equal 5. 

If one restricts to the abelian part, one has various possibilities for expressing the the- 
ory by using arguments as in [28] . The dual action, up to gauge-fixing and supersymmetric 
terms, should take the form of a very simple 5-dimensional topological term: 



J5 

Such a Chern-Simons type action is interesting, since it can be gauge-fixed in the frame- 
work of the self-duality for 2-form gauge fields [29]. Thus, the question also arises whether 




(8.3) 
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a reliable link can be established between the TQFT generated by (8.3) and the physical 
Yang-Mills theory! 

Let us finally mention the following observation. Suppose we start from an action 
with two abelian two-forms in six dimensions, /g (iS2 A dB2. The field content that one 
obtains by dimensional reduction in four dimensions will be made of 2 two-forms, 4 one- 
forms and 2 zero-forms. But the 2 two-forms are equivalent by duality to 2 zero-form in 
four dimensions. So, we obtain a field content equivalent to 4 one-forms and 4 zero-forms. 
It is striking that this is the bosonic content of the Weinberg-Salam model, with its four 
gauge fields and four scalars, which are the components of a complex SU (2) doublet. 
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Appendix A. Derivation of lattice Fokker-Planck equation 

We shall demonstrate that the discrete Markov process defined in sec. 6 is described 
by the Fokker-Plack equation of sec. 6 in the limit e — > 0. We shall derive the lattice 
Fokker-Planck hamiltonian Hpp from the transition matrix T{Ut+i,Ut), eq. (6.1) , in 
the same way that the Kogut-Susskind hamiltonian is derived from the transfer matrix of 
lattice gauge theory. We follow here the discussion of Creutz [30] . We introduce a basis of 
states labelled by group elements \U >, where U = {Ux,^}, with operators 

Ux,^^\U >= U,,^\U > . (A.l) 
We must find an operator T with the property 

<U'\T\U >=T{U',U). (A.2) 
To order e it will have the form T = 1 — eHpp, which will allow us to identify Hpp. 
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From eqs. (6.1), (6.3), (6.20) and (6.37), we have 

r(C/', C/) = AAexp ( - (e^o')"' E " ^'x^x.^t^..^)) , (A.3) 

Here H^^^ e SU{N) is given by 

Hx,i^ = exp(-eva;) exp(e/^,^) exp(e Ua;,i^fx+fi,i^U-l) , (A.4) 
and the gauge-fixed value of Ax,5 appears, namely Ax^5 = evx where 

vx = -GxT . (A.5) 

We write 

Hx,,j, = exp{ekx,^) , (A.6) 
where kx,fj, may be interpreted as the total drift force. To order e, it is given by 

kx,fi = fx,fi + {Di^v)x , (A.7) 

where is the lattice ganger covariant difference, eq. (6.30). It is the sum of the drift 
force fx,i^ = -Ex,i^S and a "gauge-fixing force" {D,^v). 

Let To be the operator whose matrix elements give the transition probability with 
total drift force kx^^ = 0, 



< U'\To\U >= jVexp ( - {eg^^)-'Y.^e Tr{I - C/',>.,^)). 



(A.8) 



We have 



< U'\T\U >=< U'\To\ exp(e/c)t/ >, (A.9) 
where exp(e/c)C/ = {exp{ekx,ij.)Ux,n} ■ We insert a complete set of states. 



< U'\T\U >= J dU" < U'\To\U" >< U"\ ex.p{ek)U > 



(A.IO) 



To order e we have 

< U"\ exp(ek)U > = S(U",exp(ek)U) = S(exp(-ek)U" , U) 

= {l-eJ2E\,,kxAU))5{U",U), (^-11) 
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where E"x,fj, is the color-electric field derivative operator defined in sec. 6, that acts on 
the variable U". Let £^a;^^ be the corresponding quantum mechanical color-electric field 
operator. In terms of operators we have shown that to order e 

T = To {l-eJ2K„KAU)) ■ (A.12) 
The evaluation of Tq may be found in [30] , with the result 

T=l + e (^o4% - KAU)) } . (A.13) 

x,m 

Thus with T = 1 — eHpp and k^^^^ = f^^^ + (I^^f)^;, and fx^/j, — —E^^f^S, we obtain the 
Fokker-Planck hamiltonian 

Hfp = - ^ -E^,^ {9oEx,fj, - {D^v)x + Ex^i^S^ , (A.14) 

x,m 

where Vx = —GxT. 



Appendix B. Equivalence of topological lattice action and lattice Fokker-Planck 
equation 

We wish to show that the 5-dimensional partition function Z, eqs. (7.46) and (7.48), 
that was obtained from the topological lattice action by integrating out all the auxiliary 
fields, is a path integral representation of the solution to the Fokker-Plank equation, d^P = 
—HppP, with Fokker-Planck hamiltonian Hpp given in eq. (A.14). 

For this purpose consider the integral defined for arbitrary 0(U') by 

0{U)= J dU'0{U')lC{U',U) 

= J dU' 0{U') ATexp ( - (4e)-i $^(5.,^ - ^Ki^f) , 



(B.l) 



X,IJ. 



where kx,^ — kx,fj,{U) is defined in eq. (A. 7), and Bx,fj, is defined by exp(ta-B"^^) — 
U' x,fj.U~jj^. It will be sufficient to show that in the limit e — this expression corresponds 
on the one hand to a discretization of the lattice action (7.48), and on the other hand 
yields 

0(t/) = (l-e4p)0(C/) . (B.2) 
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We change variable of integration by translation on the group U' = VU, with dU' = dV. 
With Vx,ij, parametrized by = exp(t(j5^^^), we have 

0{U)=Af I dBp{B)0(^eMB)u) exp ( - (4e)-i - eA;,,^)^) , (b.3) 

where p{B) is Haar measure. We again change variable of integration by Bx,^ = e^Cx,^ + 
^kx,n-, and obtain 

0{U)= j dC p{e^C + ek)0{eM^^C + ek)U)Ne^v{-\Y.Cl^) ■ (B.4) 

X,IJ, 

We now expand to order e, 



Px,ni^''C + eA;) = 1 + const x eC^^^ 



0(exp(e^C + ek)u) = (l + (e^Q^ + efc« + ^eQ^Cj,,^" ^i?J,,)o(C/), 



(B.5) 



where a sum over repeated indices is understood. When this is substituted into (B.4), the 
contribution from the expansion of p is cancelled by the normalization constant J\f, and 
one obtains eq. (B.2) , with 

^FP = -Yl {-^^'1^ + ^»,m(^))^»,m ' (B.6) 

as desired. 

We next show that the kernel in (B.l) corresponds to a discretization of the path inte- 
gral with action (7.48). In the preceding calculation the drift force k = k{U) depends on U 
but not on U', whereas the path integral (7.46) and (7.48) implicitly uses the symmetrized 
drift force 

= ^(^k{U) + k{U')) . (B.7) 

So to compare with the path integral, we express the integral (B.l) in terms of kg plus a 
correction which will turn out to be the well-known Ito term. We have 

k{U) = ks + ^(k{U)-k{U')) , (B.8) 

so, with U' — exp{B)U, we may write to the order of interest 

k{U) = kg-^YB^^^E^^^k , (B.9) 
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because the last term gives a contribution of order e as we shall show. This gives, to the 
order of interest, 

exp{-(46)-^5](5,,^-6/c,,^(C/))'}=(l-^ Yl <m^.%^.%<m) 

X exp - (4e)-^ Xll^x,^ - ^ks,x,^,f^ ■ 

Upon substituting this expression into (B.l), one obtains a Gaussian integral and Gaussian 
expectation value 

< Bl^Bl, >= 2e5,,y5^,J'''' , (B.ll) 

which is indeed of order e, as asserted. Consequently in the exponent we may replace 
B^^^By j^ by its expectation value, and to the order of interest we have 

exp{-(46)-i ^ [b,,^ - ek,,^{U)) } ^ exp{- ^ ((46)-^(S,,^ - eks,,,^f 

x,tx x,tx (B.12) 

the last term being the Ito term. Observe that with U' — > Ut+i and U — > Ut-, and 
U'x,ij,U~j^ = exp{Bj;^n), we have B^^n ~ ^Ux,nU~^^, so the last expression may be written 
formally as 

exp{-e J] {^{Ux,,JJxi, - ^Kx,^? + 2^S,m^S,m)} • (B.13) 

X,iX 

We have 

K,,K,, = - GIGIT , (B.14) 

and the last expression corresponds to the action (7.48), as asserted. 
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